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Abstract 

The Bfc parameter is computed in quenched lattice QCD with Wilson twisted mass fermions. 
Two variants of tmQCD are used; in both of them the relevant AS" = 2 four-fermion 
operator is renormalised multiplicatively. The renormalisation adopted is non-perturbative, 
with a Schrodinger functional renormalisation condition. Renormalisation group running is 
also non-perturbative, up to very high energy scales. In one of the two tmQCD frameworks 
the computations have been performed at the physical i^T-meson mass, thus eliminating the 
need of mass extrapolations. Simulations have been performed at several lattice spacings and 
the continuum limit was reached by combining results from both tmQCD regularisations. 
Finite volume effects have been partially checked and turned out to be small. Exploratory 
studies have also been performed with non-degenerate valence flavours. The final result for 
the RGI bag parameter, with all sources of uncertainty (except quenching) under control, is 
Bk = 0.789 ± 0.046. 



1 Introduction 



Indirect CP-violation in K ^ tttt decays is expressed in terms of the e-parameter. Its 
known experimental value, combined with theoretical input from neutral ii'-meson 
oscillations, defines a hyperbola in the complex plane of the unitarity triangle (for 
a review see ref. [1]). The theoretical prediction of the oscillation amplitude is 
obtained in an operator product expansion framework, as the product of a single 
Wilson coefficient (known to NLO in perturbation theory [2,3]) and the matrix 
element (K°|0^'^=^|K°) of a dimension-6, AS = 2 operator 

0^^=2 = [,-^^(l-^5)d] [s-7^(l-75)d] 

= Ovv+AA - OvA+AV • (1-1) 

Here s and d are strange and down quark fields. Within the square brackets spin 
and colour indices are saturated. In obvious notation, the operators Oyv+AA and 
OvA+AV are the parity-even and -odd parts of O^'^^^. For historical and techni- 
cal reasons, this matrix clement is usually normalised by its value in the vacuum 
saturation approximation (VSA). It is thus expressed in terms of the Bk parameter: 

_ (A-"|Ovv+aa|A-") (A'"KA-v+aa|/<°) 2^ 
^ - (KO|Ovv+AA|i^o)y5A 8 ' ^ - ) 

with rriK the X-meson mass and Fk its decay constant. The above expression 
involves the matrix element of the parity-even operator Oyv+AA only. The parity- 
odd matrix element {K^\0\x+KS/\K^)vSA vanishes identically in QCD. 

Bk is an intrinsically non-perturbative quantity, which can be computed in the 
lattice rcgularisation of QCD. The results of the various recent Bk measurements 
have been summarised in ref. [4] (see also references in this review). Besides quench- 
ing, which is a still uncontrolled source of systematic error, the most important 
source of uncertainty in these computations arises from the operator renormalisation. 
In schemes which respect chiral symmetry, the operator Oyv+AA is multiplicatively 
renormalisable. Thus, the corresponding physical matrix element can in principle 
be accurately obtained by a series of non-perturbative (lattice) measurements of the 
bare matrix element and the operator renormalisation constant, at several values of 
the lattice spacing, followed by a continuum limit extrapolation. This is the case 
of Ginsparg- Wilson fermions, where chiral symmetry-breaking effects arc negligible. 
With Wilson fermions, the breaking of chiral symmetry induces mixing with four 
other dimension-6 operators [5-9] in the parity-even sector: 

4 

(OR)yv+AA = -^vv+AA(5o,aAi) Ow+aa(5o) + Xl^»(fo)^'(fo) ■ (^■^) 

1=1 
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The operators Oi{go) belong to different chiral representations than Oyv+AA- The 
mixing coefficients Aj(^o) are finite functions of tlic bare coupling, while the renor- 
malisation constant ,^vv+Aa(5'0 7 fJ/^) diverges logarithmically. This mixing pattern 
is the reason behind the relatively poor precision of Bk, measured with Wilson 
fermions. 

Two proposals have attempted to resolve this issue. They are both based on 
the observation [7, 9] that, even in the absence of chiral symmetry, the parity-odd 
operator OvA+AV is protected from mixing with others, by discrete symmetries: 

(Or)va+av = -^VA+Av(5o, a/u)OvA+Av(5o) ■ (1-4) 

The first proposal [10] consists in obtaining the physical — K'^ matrix element of 
Ovv+AA from a correlation function of the renormaliscd operator OvA+AV? through 
axial Ward identities. The method has been put to test in rcf. [11]. The Bk estimate 
of this method turned out to be compatible with the result of the standard com- 
putation (with operator subtractions). Unfortunately, the correlation function of 
Ova+AVj being the product of four composite operators, turned out to be noisy and 
the error was as large as the one of the standard computation. Thus, this method 
is successful in eliminating an important source of systematic errors (operator sub- 
traction) at the cost of increased statistical fluctuations. 

In this work we implement the second proposal [12], which is based on the 
addition of a so-called "twisted mass term" to the Wilson fcrmion action. This entails 
loss of parity and partial loss of flavour symmetry at finite lattice spacing, recovered 
in the continuum limit. On the other hand, some renormalisation properties are 
greatly ameliorated in the twisted mass QCD (tmQCD) formalism. The relevant case 
for us is the renormaliscd (X*'|Ovv+AA|-f^^) matrix element, which in the tmQCD 
formalism may be expressed in terms of the parity-odd operator OvA+AV- The 
tmQCD action differs from the standard one by an additional soft term, which does 
not modify renormalisation properties in mass independent renormalisation schemes. 
In particular, the operator OyA+AV remains multiplicatively renormalisable, with the 
same renormalisation constant and running as with Wilson fermions. Thus finite 
subtractions are avoided in the tmQCD determination of Bk- 

With respect to previous studies, we have introduced several new features: 

• The bare matrix elements have been obtained in two distinct lattice tmQCD 
formulations. In the first formulation, the down quark is twisted, with a twist 
angle a = tt/2, while the strange quark is a standard (untwisted) Wilson 
fermion. In the second formulation we have a twisted flavour doublet of down 
and strange flavours, with twist angle a = 7r/4. Having two independent Bk 
estimates allows a better control of the continuum limit extrapolations. We 
have implemented Schrodinger functional (SF) boundary conditions for these 
computations. 
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• Quenched simulations with standard Wilson fermions are carried out at heavier 
quark masses, due to the presence of exceptional configurations. Since our 
tmQCD variant with a = tt/A is free of exceptional configurations, Bk can 
be computed at the physical if -meson mass value. Thus the error due to 
extrapolations from higher masses is eliminated. 

• In ref. [13] the operator OvA+AV has been renormalised non-perturbatively 
in several SF renormalisation schemes and its Renormalisation Group (RG) 
running has been computed non-perturbatively from low energies up to scales 
of several tens of GeV. We have used these results in order to obtain the 
Renormalisation Group Invariant (RGI) counterpart of Bk- Our result is 
essentially free of any uncertainty due to Perturbation Theory (PT). 

• Although our main results have been obtained with degenerate down and 
strange quarks, we have also performed some exploratory studies with non- 
degenerate flavours. Also, some simulations have been performed in order to 
probe finite volume effects. 

Our result confirms earlier ones, obtained from other lattice regularisations (domain 
wall, overlap, staggered) albeit with a major control of the various sources of error. 
An interesting discrepancy with the result of the most recent detailed study of Bk 
with Wilson fermions (see ref. [11]) is also resolved. 

The paper is organized as follows: In sect. |21 we review the tmQCD basics, 
emphasising the cases of interest. We show how, by choosing either a twist angle 
a = 7r/2ora = 7r/4, we can map the operator Oyv+AA onto OvA+AV) thus achieving 
multiplicative renormalisation without subtractions. In sect. 01 we review how four- 
fermion operator matrix elements can be obtained from SF correlation functions. In 
sects. Eland El we present our raw results for tmQCD with a = 7r/2 and a = 7r/4 
respectively. In sect. El we compute the RGI- extrapolated to the continuum 
limit. Several more technical aspects are discussed in the appendices. 

Preliminary results of our work have appeared in refs. [14, 15]. Based on a 
similar approach, the B^ — B^ matrix element, with a static bottom flavour, is 
currently under study [16,17]. 

2 Twisted mass QCD and weak matrix elements 

Twisted mass QCD has been designed to eliminate exceptional configurations in 
(partially) quenched lattice simulations with light Wilson quarks [12]. In its original 
formulation, it describes a mass-degenerate isospin doublet if) of Wilson quarks for 
which, besides the standard mass term, a so-called twisted mass term ijiq'ip'y^T'^ip is 
introduced. The properties of tmQCD have been studied in detail in [12], where, in 
particular, its equivalence to standard two-flavour QCD has been established. We 
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discuss here two extensions of this framework, which are suited to the extraction of 
Bk. 

2.1 Twisted mass QCD with twist angle 7r/2 

The first variant has aheady been mentioned in ref. [12], and corresponds to a 
Euchdean continuum action of the form 

Sf = j d^x [i}{x){Ip + mi + iiia5T^)ipix) + s(x)(^ + m,)s(x)] . (2.1) 

Here, the two Ught flavours form an isospin doublet = (u, d), is a Paiili matrix 
and fii and m/ arc the twisted and standard quark mass parameters. For the hght 
quark sector the physical quark mass is given 

Ml = .f^q + ?i , (2.2) 
and the twist angle is defined by 

tana = -^. (2.3) 

mi 

The considerable advantages of the choice a = 7r/2, for which = and Mi = /ii, 
will be discussed below. We will refer to this case as "fully twisted", since the 
physical quark mass M; is determined by the twisted mass parameter alone. 

The equivalence of this theory to standard QCD, established in [12], is based 
on axial transformations of the quark fields and corresponding spurionic transforma- 
tions of the mass parameters. Let us relabel the fields and masses of the above theory 
by il)'^'^\mf^ and ^if^. The corresponding quantities in standard QCD (which is 
tmQCD at q; = 0) are given by ij;^^^ , mf^ and = 0. The two theories are then 
related by the axial field transformations 

^(a) ^ ^(0) ^ ^(") ^ ^(0) = ^(")i?(a) , (2.4) 
and the transformation of the mass parameter 

rnf' = m["^ cos(a) + sin(a) , (2.5) 

where 

i?(a) = exp{i75ar3} . (2.6) 
In tmQCD we denote Euclidean correlation functions by 

<0[V', V^]>. . = Z-^ [ 0[^,^]e-^, (2.7) 

^ ' ^fields 
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where Olipjip] denotes some multilocal gauge invariant field. The fermionic part 
of the action is given in ()2.1|) . The relation between standard QCD and tmQCD 
correlation functions, are then expressed as 

(0[V',?]>(,) = (O[i?(-a)^,?i?(-a)]>(o), (2.8) 
(0[V',?]>(o) = (0[i?(a)V',W«)]>(„). (2.9) 

Hence, a given correlation function in tmQCD with twist angle a is interpreted as 
the linear combination on the r.h.s. of Eq. ()2.8|) . If instead we are given a standard 
QCD correlation function, then Eq. (|2.9|) tells us how it is represented in tmQCD 
at twist angle a. 

In the formal continuum framework, all these relations between tmQCD and 
standard QCD quantities are readily obtained from the chiral transformations (|2.4j) . 
Their extension to the lattice regularised theory with Wilson fermions, which is the 
lattice regularisation of choice in the present work, is more intricate. As shown 
in [12], the same relations are realised between the renormalised quantum field theo- 
ries, provided the renormalisation procedure is set up with some care. Moreover, the 
definitions (|2.2|) and ()2.3() are understood to hold for suitably renormalised masses 
Mr^i, mji^i and /xr,«, such that the twist angle a is free of renormalisation. 

2.1.1 Relations between composite fields 

The axial transformations ()2.4() induce a mapping between composite fields. For the 
quark bilinear operators 

Sij = ipiipj, Pij = i>a5il^j, A^^ij = tpa^.-y^ipj, v^^ij = ■ipa^.ipj, (2.10) 

we have 

(2.11) 
(2.12) 
(2.13) 
(2.14) 

For the four-quark operators under consideration we have: 

^VV+AA = COS (a) Oyv+AA - ^ sin (a) OyA+AV ' (^-l^) 
'^VA+AV = cos(a)0^'^^Ay - isin(a)0j'^_^AA. (2.16) 



o(0) 


= COS 


{2) '^sd 




p(0) 
^sd 


= cos 


(a\ p(a) 
2 ; ^sd 


-zsin(f)5S\ 


AO) 

^fj.,sd 


= cos 


(a\ 4(") 
{2) ^fi,sd 




fi,sa 


= cos 


{2) ^fj.,sd 
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The case of interest is q = 7r/2, for which we have 

•^VV+AA ~ "^VA+AV "-^VA+AV ~ *'-^VV+AA • l^^-iOJ 

The above formal expressions imply the following relation between renormalised 
operator matrix elements: 

{K'\ (OR)^°i+AA \K') = -i{K^\ (OR)tr/i)AV ■ (2.19) 

The advantage of computing Bk from the matrix element on the rhs (i.e. by per- 
forming a tmQCD computation), is that OyA+AV is multiplicatively renormalisable, 
whereas the matrix element on the Ihs involves Oyv+AAi which requires additive 
renormalisation. 



2.2 Twisted mass QCD with twist angle 7r/4 

We also implemented a second variant of tmQCD, in which the down and strange 
quarks are grouped into a flavour doublet = {s,d): 

Sf = j d^x ['ilj{x){]p + mi + i-f5r'tii)Hx)+uix)ilp + rnu)uix)] . (2.20) 

This formulation of the theory is suitable for Bk studies with degenerate light 
and strange valence flavours. In the quenched approximation, the first term alone 
is adequate for Bk simulations, since only strange and down valence quarks are 
involved. In this approximation it is straightforward to also accommodate non- 
degenerate valence flavours, by introducing diagonal 2x2 mass matrices in place of 
mi and /i/ (see ref. [18]) for details^). 

The general relationships of subsect. 12. H defining total quark mass, twist angle 
and chiral rotations between fermion fields, mass parameters and correlation func- 
tions, hold also in this case. Once more, they are valid both formally and between 
renormalised quantities. 

^In the present work, simulations with non-degenerate valence quarks have only been carried 
out with twist angle a = tt/2. 
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o(0) 


_ qia) 
— '^sd ' 




p(0) 


_ p(") 

~ ^sd ' 




H,sd 


= cos(a)^2<i 


-zsin (a)V;J°)i, 


y(0) 
H,sd 


= cos(a)V;J°); 


-zsin (a) ^J;^]^. 



^.S.i Relations between composite fields 

In the present tmQCD formulation, the chiral rotations between composite fields 
become 

(2.21) 
(2.22) 
(2.23) 
(2.24) 

For the four-quark operators of interest we have: 

0t°i+AA = cosi2a)0^^i^j,^-ism{2a)0^4^^^, (2.25) 

OvA+AV = cos(2a)o(,^+AV-^sin(2a)o(r^+AA- (2-26) 
At a = 7r/4 we have 

^U = %^if-'Vi:Jfi. cm 

Wv+AA — ^Wa+AV ' '^VA+AV ~ ^Wv+AA ■ l^z.zo; 

The formal above expressions imply for the renormalised WME of interest 

(Or)Vv+aa = -^{K'\ {Or)'^IX • (2-29) 

Here, as in the tt/2 case, we see that the QCD four-fermion WME of interest is 
mapped onto a tmQCD WME which involves the multiplicatively renormalisable 
four-fermion operator OvA+AV- So from the renormalisation point of view, the 7r/4 
version is equally advantageous. 



2.3 Flavour symmetry, strangeness and mass degeneracy 

The two tmQCD formalisms exposed above have distinct characteristics, which merit 
some discussion. The it/2 case refers to a light flavour doublet, while the strange 
quark is regularised in the standard way. Thus simulations may be naturally per- 
formed with non-degenerate down and strange flavours. In particular, simulations 
may get close to the physical situation, since the twisted light quarks are protected 
from exceptional configurations, whereas the non-twisted strange quark is heavy 
enough to remain unaffected by this problem. 

At this point we note that most quenched computations of Bk have remained 
with mass degenerate down and strange quarks, simply because quenched chiral 
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perturbation theory indicates the appearance of potentially dangerous quenched 
chiral logs as soon as one departs from this situation [19-22]. In this unphysical 
situation, the problem of exceptional configurations re-appears for the strange quark. 
Thus one is forced to stay with fairly massive pseudoscalar mesons of about 600 MeV, 
where the problem of exceptional configurations is believed to be absent, at least 
with a non-perturbatively 0(a) improved action [23]. In quenched computations 
with degenerate strange and down masses, one must therefore compute Bk with the 
if —meson tuned to several values above its physical mass and then extrapolate to 
the physical point. 

While in the a = 7r/2 scenario some tuning of the quark mass parameters is 
necessary to achieve mass degeneracy between the down and the strange quark (see 
below), this unphysical situation is naturally obtained in the 7r/4 case where down 
and strange quarks form a flavour doublet. The problem of exceptional configura- 
tions is now absent. This enables the computation of Bk with degenerate valence 
quarks tuned so as to have the if-meson at its physical mass value, avoiding extrap- 
olations from heavier masses. 

More generally, it must be realised that partial loss of flavour symmetry at finite 
lattice spacing is the price one has to pay for the attractive features of tmQCD. Sur- 
prisingly large flavour breaking effects have been observed with maximally twisted 
Wilson quarks at Cgw = in refs. [24, 25] . We have also looked at flavour breaking 
effects and found that these are reasonably small and rapidly decreasing towards 
the continuum limit (cf. sect.|lj). The question whether these different findings are 
due to differences of the lattice action, parameters or the choice of observables will 
be the subject of further study. 

2.4 Improvement considerations 

Lattice quantities are computed at several fixed UV cutoffs and, after they are 
renormalised, they are extrapolated to the continuum limit. Symanzik improvement, 
if applied, ensures a better control of the extrapolating procedure. It involves adding 
0(a) counterterms both in the lattice action and the operators. In this work we will 
be using the Clover improved action [26] and improved currents (in the spirit of 
ref. [27,28]). 

Use of the Clover improved action implies that our meson mass estimates are 
subject to finite spacing effects which are O(a^). In order to remove 0(a) cutoff 
effects from Bk, we would also need to improve the relevant matrix elements of 
the four-fermion operator and the axial current in the tmQCD regularisation. The 
improvement pattern of dimension-3 quark bilinear operators in the quenched ap- 
proximation is given in Appendix El For tmQCD with untwisted (strange) and 
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twisted (down) quarks, these results give for the currents 

{AR)^^sd = Za[1 + lbj^amg^s][A^^sd + CAad^Psd - i^afiibj^V^^sd] , (2.30) 
{VR)f,,sd = ZY[l + ^b^amg^s][V^,sd + '^yo,duTf,^,sd-ila^iKAf,^sd], (2.31) 

where the subscript R refers to renormaUsed quantities and d stands for the sym- 
metric lattice derivative. These renormalised currents combine as in Eq. I|2.17|) to 
give the improved axial current in the a = tt/2 case. For tmQCD with the down 
and strange quarks in the same (twisted) doublet, the results of [29] can be directly 
taken over (with the up quark replaced by the strange one). We thus have 

(ARj^^sd = Za[1 + bj,amq\[Af,^sd + CAadfj,Psd - iafiibj^V^^sd] , (2.32) 

{VR)^^sd = Zy[1 + ftvomq] [V^^sd + cyaduT^uM - iamKA^^sd] ■ (2.33) 

These renormalised currents combine as in Eq. (|2.27|) to give the improved axial 
current in the a = 7r/4 case. The renormalisation constants and improvement 
coefficients used here may be found in Appendix 1X1 

The improvement of four-fermion operators is a far more difficult procedure. 
Though feasible in principle, it is rendered impractical by the proliferation of coun- 
terterms. We will hence not proceed in this direction.^ 

Following these considerations, we have always used the Clover action in our 
simulations. The four-fermion operator is left unimproved. The implications of 
current improvement to the 0(a) effects of Bk require further discussion. We will 
repeatedly revert to this issue at later stages. 

2.5 Tuning of quark masses 

With the continuum actions ()2.1|) and (12.20(1 the twist angle is directly determined 
by the ratio of the mass parameters mi, m in the action. This is no longer the case 
once tmQCD has been regularised with Wilson type quarks. Rather, the tuning of 
the twist angle a to the preferred values tt/2 and 7r/4 requires the implementation 
of Eq. (|2.3() for renormalised mass parameters. We discuss this issue in some detail. 

We denote the subtracted bare quark mass for Wilson fermions by omq = 
1/(2k) — l/(2Kcr), K being the hopping parameter. Whenever we need to identify 
the quark flavour /, we will denote the corresponding quantities by am^j and kj. 
For the untwisted strange quark in the 7r/2 formulation, the renormalised quenched 
quark mass is given by 

niR^s = Z^lrriq^sil + braarriq^s)] , (2.34) 
^The proposal of ref. [30] is an interesting alternative, which will also not be pursued here. 
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while the Ught quark masses renormaUse as follows: 

rriR^l = Zm[m.q,/(1 + bmarriq^i) + bmafif] , 

^J■R,l = Zf,fii{l + b^amq^i) . (2.35) 

The above expressions are valid up to O(a^) corrections. For the strange flavour, 
Eq. (|2.34)) is the standard relationship between renormalised quark mass and Ks- For 
the light sector, fixing the twist angle to it/2 + 0{a?) amounts to setting m-^^i = 
with ^Yi^i positive. In other words, once a value a/x; is chosen, we tune so as to 
satisfy 

arriq,; = -6m(aw)^ . (2.36) 

The 7r/4 case is somewhat less trivial. In order to ensure that the twist angle 
acquires this value to O(a^), we must impose /Ur = mR to this order in the cutoff 
effects (flavour indices are suppressed, as the formalism applies to degenerate down 
and strange quarks). In terms of eqs. ()2.35|) this means 

amq = — ^a/xq{l + [-r^ip^ - h^) - ZZAbra]aHq} , (2.37) 

with Z = Z^/{Zf^Z\). For a given choice of afiq and tuned so that arua 

satisfies the above relation. 

In Appendix^ we collect the known results for the renormalisation constants 
and improvement coefficients required above. 

3 SF Correlation functions at large time separations 

We now derive explicit expressions for the representation of Schrodinger functional 
correlation functions in terms of intermediate physical states. For quark bilinear 
dimension-3 operators (e.g. the pseudoscalar density, or the axial vector current) 
this has been discussed in ref. [31]. Here we recapitulate the derivation and generalise 
it, in a straightforward manner, to the four-quark dimension-6 operator of interest. 
We then discuss how this formalism extends to the case of tmQCD. 

3.1 Quantum mechanical representation of the Schrodinger functional 

The QCD Schrodinger functional is defined as the QCD partition function in Eu- 
clidean space-time, with quark and gluon fields obeying periodic boundary condi- 
tions in space (with period L) and Dirichlet boundary conditions in time at the 
hypersurfaces xq = and xq = T. We assume homogeneous boundary conditions, 
i.e. the spatial components of the gauge potentials and the Dirichlet components of 
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the quark and anti-quark fields are taken to vanish at the two time boundaries. The 
Schrodinger functional can then be written as [32,33] 

Z = (io|e-™P|io), (3.1) 

where |io) is a state which is implicitly defined by the Dirichlet conditions, and 
carries the quantum numbers of the vacuum state. The presence of P implies a 
projection onto the subspace of gauge invariant states, and H denotes the Hamilton 
operator of QCD formulated on a torus of volume . It is the existence of the 
Hamiltonian operator which allows for the definition of a time-zero quantum me- 
chanical Hilbert space and the corresponding operator representation of Euclidean 
correlation functions. On the lattice, aH is defined as the negative logarithm of the 
transfer matrix and it is Hcrmitian provided the transfer matrix itself is Hermitian 
and positive. This is indeed the case with the standard Wilson gauge action and 
unimproved Wilson quarks with both standard and twisted mass terms [29,34]. In 
principle, this property is lost in the presence of 0(a) improvement terms in the 
action. However, the ensuing unitarity violations are usually considered harmless 
since they occur close to the cutoff scale and thus do not affect the physics at low 
energies [35]. 

The quantum mechanical representation of correlation functions then proceeds 
via the introduction of a set of (gauge invariant) eigenstates of H, 

|n,(?), n = 0,l,... , (3.2) 
M|n,g)=4«)|n,g), (3.3) 

with normalisation {n',q'\n,q) = 5n^n'^q,q'- Here, g is a multi-index comprising all 
quantum numbers corresponding to the exact lattice symmetries, and n enumerates 
the energy levels in the channel specified by q. We do not indicate the momentum 
of the states |n, g), since we will always be working with correlation functions that 
project states with vanishing (spatial) momentum. 

Note that in general the set of conserved lattice quantum numbers q is smaller 
than in the continuum, due to the explicit breaking of symmetries by the lattice 
regularisation. Standard Wilson quarks have the nice property to conserve an exact 
SU(A^f) flavour symmetry, besides parity and charge conjugation symmetries. The 
only difference to the continuum classification of particle states then consists in 
the breaking of rotational symmetries by the lattice regularisation, which implies 
that angular momentum is not a good quantum number in general. However, this 
mainly affects higher spin states, and is irrelevant for the pseudoscalar meson states 
at hand. In the following wc first give an account of this continuum like situation 
with standard Wilson quarks before turning to the more complicated case of tmQCD. 
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3.2 Specific cases of SF correlation functions 

For the calculation of Bk, we are interested in specific correlation functions of gauge 
invariant dimension-3 quark bilinear operators X (where X may denote scalar and 
pseudoscalar densities S, P or currents V^, A/^), the four-fermion operator Ovl"^ 
and the dimensionless boundary quark fields 

^ds = ^Yl Cd(y)75C.(z) , = ^ E Cd(y)75CU^) , (3.4) 

which have been discussed in [28] . In terms of these quantities we define the gauge 
invariant correlation functions 

^3 

hi^o) = -Yi^sd{x)Ods), (3.5) 

where (...) denotes the usual path integral average. Upon specifying X = S, P,Vfj,, 
we denote the above correlation as /x = fs, fp, fv, Ia respectively. Similarly, for a 
AS" = 2 four-fermion operator Y we define the correlation function 

Fy(xo) = L'^iO'a, Y{x) Ods) . (3.6) 

Upon specifying Y = Ova,Ow)Oaa we denote the above correlation as = 
Fva,FvvjPaa respectively. Finally, in order to obtain properly normalised hadronic 
states one needs to consider the boundary-to-boundary correlation function 

/i = -l{Oi,Ods). (3.7) 

Clearly, in the standard lattice formulation the parity-odd correlation functions 
/s,/v and FyA vanish identically, but this will not be the case in tmQCD, where 
the corresponding fields will be re-interpreted according to the discussion in sect. 2. 
The correlation functions /x have the quantum mechanical representation [31]. 

/x(xo) = ^-^y (io|e-(^-"°)=Xe-^»=P|iK), a < xq < T - a , (3.8) 

where X is the corresponding time-independent operator, and the state jix) has 
the quantTim numbers of the K-meson with momentum zero. Analogously for the 
correlation function Fy = i^vv^^AA we obtain 

Fy(xo) =.Z-^L6(iK|e-(^-^°)^Ye-^°^P|iK), a<xo<T-a, (3.9) 

i.e. the operators O^s and O'^^ carry the quantum numbers of a K'^ by construction. 
We also have 

/i=Z-ii(iK|e-™P|iK). (3.10) 
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The asymptotic behavior of fx{xo) (with X = Aq,P), Fy(xo) (with Y = Ovv+aa) 
and fi for large values of both xq and T — xq (with L unspecified at this stage) is 
as follows 

/x(xo) « f p(0,0|X|0,i^)e-^o'"^ X {l + 77fe-^«^ + r?V^^""°^™°}, 
Fy(xo) ~ L6|pp(0,K|Y|0,K)e-'^""^, 

/i « ip^e-^--, (3.11) 
where we have introduced the ratios 

(0,i^|iK) 



(0,0|io) ' 



(3.12) 



^ _ (0,0|X|l,i^)(l,K|iK) 

" (0,0|X|0,i^)(0,K|iK)' ^ ' 

^ (io|l,0)(l,0|X|0,i^) 

(io|o,o)(o,o|x|o,i^) ■ ^ ' 

The energy difference uiq = Ei'^ - is the mass of the 0'*'"'' glueball and A = 
e[^^ — Eq^^ is an abbreviation for the gap in the ii'-meson channel. We have 
dropped contributions of higher excited states which decay even faster as xq and 
T — Xq become large. 

Considering the special case of /a, one finds that it is proportional to the matrix 
element (0, 0| Ao|0, -fC), which is related to the kaon decay constant Fk through 

Za(0, 0| Ao|0, K) = FKmK{2mKL^)~^/'^ . (3.15) 

Here, Zx is the renormalisation constant of the isovector axial current, and the 
factor (2mj<--^^^)~^^^ takes account of the conventional normalisation of one-particle 
states^. In our convention the experimental value of the pion decay constant is 
132 MeV. 

Eq. ()3.11|) is used to determine mj^, while the Kaon decay constant Fk may be 
conveniently extracted from the ratio 

ZAfA{xo)/VTi « iF^(mKL3)V2e-(^»-^/2W 



X {l + r/f e-^°^ + r/V^^""'"^""'} • (3-16) 



'We denote conventionally normalised one-Kaon states by \K). 
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Similarly, the AS" = 2 matrix element of the operator Y can be determined from the 
ratio 

fl^.L«(0,A-|Y|0,A-)=^(A1Y|«,, 

where in the last equation the factor {2mKL^) again refers to the conventional 
normalisation of one-particle states \K) and \K) in infinite volume. Finally, the 
bare Bk can be extracted from ratios of the form 

^y(^o) _ (0,K|Y|0,K) ^^^^^ 



;[2/A(xo)][2/;(r - xo)] 1(0, 0| Ao|0,i^)(0, K\ Ao|0, 0) 

{K\Y\K) 



Zi 



3^K'"-K 



where we define 



'3 



f^(T - xo) = -^{0',,Xasix)) , (3.18) 

with an analogous asymptotic expansion to that of /x (note that, for this correlation, 
the dominant exponential decay is exp[— (T — xo)mK])- The above formulae show 
explicitly how masses and matrix elements can be obtained from Schrodinger func- 
tional correlation functions. A discussion on the practical advantages of this method 
of extraction of hadronic masses and matrix elements can be found in ref. [31]. 



3.3 SF correlation functions and tmQCD 

Twisted mass QCD requires two modifications to the above framework. First, the 
fields and the symmetries must be re-interpreted through the axial field transforma- 
tion as discussed in section 2. Unfortunately a direct comparison of SF correlation 
functions between tmQCD and standard QCD is not possible with our current set- 
up, i.e. the relations (|2.8() and (|2.9() . connecting tmQCD and QCD renormalised 
correlation functions in an infinite volume, do not hold for SF correlation functions 
like the ones defined in Eq. ()3.8() and Eq. ()3.9|) . For this to work out one would need 
to chirally rotate both the quark mass terms and the quark boundary projectors at 
the same time (see [36] for work in this direction) . As we do not change the fermionic 
boundary projectors, the SF correlation functions for tmQCD and standard QCD 
are inequivalent even in the continuum limit. For the quantum mechanical repre- 
sentation this implies that the initial and final states with quantum numbers of the 
vacuum or a kaon state are not the same in both cases, although we will keep the 
same notation. However, the operator relations in tmQCD and standard QCD, in- 
tended as equations between renormalised matrix elements of physical states, remain 
valid. 
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Second, the exact lattice symmetries in tmQCD, which are relevant for the 
classification of excited states in a given channel, are less restrictive, as part of the 
flavour symmetry group is broken, as well as parity. In contrast to standard Wilson 
quarks one thus has to deal with excited states which have the wrong continuum 
quantum numbers, but share all the lattice quantum numbers with the state of 
interest. A prominent example is the appearance of the neutral pion (or kaon, 
depending on the definition of the twisted doublet), which has the same lattice 
quantum number as the vacuum state. Therefore, it is always possible to have an 
excited state consisting of the hadron of interest together with a zero-momentum 
neutral pseudoscalar meson. On the other hand, it is clear that the corresponding 
matrix elements with these states are a pure lattice artefact suppressed by some 
power of the lattice spacing. This means that if one were to analyse the hadron 
spectrum after having taken the continuum limit of (ratios) of correlation functions, 
these additional states would not play any role. However, this procedure being 
somewhat impractical, any analysis at fixed lattice spacing must deal with this 
problem. 

From the above discussion and Eqs. iimi.iHll.mi'M:^.^^!^ we readily 
conclude that in tmQCD (with a = 7r/2,7r/4), Bk can be computed as the asymp- 
totic limit of the ratio 

Rb = J-q iZvA+AvFyA+AV _ 

Y[^a/a(2;o) - iZYfy{xo)][ZAfA{T - xo) - iZyf^^iT - xq)] 

The numerator of the above ratio involves the four-fermion operator which, as ex- 
plained previously, cannot be readily improved. Thus our Bk lattice estimates are 
affected by 0(a) finite cutoff effects. For the denominator of Eq. ()3.19|) we use the 
following expressions: 

ZA[fA{xo) + CAadofp{xo)] - iZyfy{xo) , (3.20) 

^a[/a(2;o) + CAadofp{xo)] - fZv/v(xo) . (3.21) 

Comparing the above with Eqs. 1)2. 3012. 3112. 3212. 33() . we note that the vector cur- 
rent counterterm proportional to ct vanishes due to invariance of /t under spatial 
translations. Moreover, we have omitted the improvement counterterms propor- 
tional to 6v,A and 6v,A- The reason is that these factors are the same for both 
currents at tree level (6a = by = 1/2 and bA = by = 0). It can then be easily shown 
that tree- level 0(a) improvement of the numerator FyA+AV would require the same 
terms as those of the denominator. Omitting these counterterms from the denom- 
inator thus amounts to a consistent cancellation of O(amq) effects at tree level in 
the ratio Rb, while there are no 0(a/iq) effects to this order of perturbation theory. 
We have nevertheless checked that, in practical simulations, the contribution from 
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the counterterms proportional to 6A,fcv)^A and by in the denominator is indeed 
neghgible. 

Concerning excited states with the wrong continuum quantum numbers we 
discuss for concreteness the asymptotic behaviour of a suitable combination of 
correlation functions, from which the decay constant Fx is extracted, ignoring 
improvement-related issues: 

^[ZAh{xo)-iZyMxo)] « |^p(0,0| A(°)|0,K)e^^°^'"- 
x {l + r?f e-^o^ + r?^e-(^-^«)'"G + r/X^g"^"^- + 7?X°e~^^"'^'')'"-« } • (3.22) 

Here the superscript in the axial current operator reminds us that the renor- 
malised operator relations corresponding to Eqs. ()2.17I2.27|) have been used. As 
these relations become exact only in the continuum limit, besides the rj^ and rj^ 
of Eqs. ()3.13|) and (|3.14|) . we have now included the lowest lying parity- violating 
contributions 

^<r^ ^ (0,0|A(°)|0,a,,)(0,a^|iK) ^3^3) 
^ (0,0|A(°)|0,K)(0,i^|iK) 

,f = (io|0^-°)(AO|Ar|0,K)^ ^3^^^^ 
(io|0,0)(0,0|A[,°>,K) 

These essentially consist of: 

• The projection of the lowest lying scalar state with a net s- and d-flavour 
content, denoted by \0,aK) from the xq = time-boundary. It has a mass 

niai^ with Act = TTT-CT — ruK ■ 

• The projection of the lowest lying pseudoscalar state with vacuum (zero) net 
flavour content, denoted by |0,7r'^) from the xq = T time-boundary. It has 
mass m^o.'^ 

In standard QCD, the relation rj'^ = rfp^ = is rigorously satisfied on the grounds 
of parity conservation. As tmQCD breaks parity, these coefficients are cutoff effects 

*Some care is required in the interpretation of this notation. Our 7r/2 tmQCD consists of a hght 
twisted quark doublet ijj = {u , d) and an untwisted strange quark. So what is meant by |0, 7r°) 
is a state generated from the vacuum by the operator tp'yzT^ip; i.e. a pion. Our 7r/4 tmQCD, 
consists of a twisted quark doublet ip = {s , d); thus we have a theory in which physical strangeness 
is described by J3 = ±1/2 isospin quark states. We therefore have three Kaons, corresponding to 
isospin values I3 = —1, 0, +1. What is meant by |0, tt") is again a state generated from the vacuum 
by the operator ip^^r'^ip, which is an Is = Kaon. 
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which only vanish in the continuum hmit. Based on our knowledge of QCD spec- 
troscopy, we expect the scalar state |0,crx) to be significantly heavier than \0,K), 
so it should have an exponentially vanishing contribution. On the other hand, the 
state |0, tt''), be it a true neutral pion in the 7r/2 theory or an J3 = Kaon in the 
7r/4 theory, is lighter than the glueball. Although the associated matrix element, 
being 0(a), is diminishing as we approach the continuum limit, the exponential de- 
cay predominates on that of the glueball, so that the net effect may be comparable 
(or even dominant) to that of the vacuum decay term. In our numerical analysis we 
obtained a rough estimate of the mass and prefactor or the first excited state, so 
that its contribution to the effective mass of the ground state could be estimated. 



4 Lattice results from tmQCD at a = n/2 

The quenched Bk parameter has been computed in tmQCD with a twist angle 
a = tt/2 at four values of the lattice spacing in the range o 0.06 — 0.09 fm, 
with the spatial directions extending from L ~ 1.4 to -L ~ 1.9 fm. Time ranges 
from = to = T, with T/L ~ 2.3 — 3.0. The three mass parameters (i.e. the 
standard hopping parameter Kg for the strange quark, the twisted hopping parameter 
Kd and the twisted mass /i^ for the down quark) are tuned as discussed in sect. HI so 
as to keep the two quarks degenerate. In order to avoid exceptional configurations 
in the untwisted strange sector, the K- meson mass is kept in the range rriK ~ 
640 - 830 MeV; for a discussion on the presence of exceptional configurations, see 
Appendix O The parameters of our runs are displayed in Table ^ Following 
ref. [37], we will express our results in units of the length scale vq ~ 0.5 fm. For the 
relationship between rg and the lattice coupling /3 (which fixes the lattice spacing 
a), see Appendix 1X1 

The effective pseudoscalar meson mass is computed from the ratio 



aM||(xo) = iln 



/AR(a^o - a) - ih^jxp - a) 
JaA^o + a) - i/vH (xo + a) 



(4.1) 



where /ar , /vr are the correlation functions constructed from the improved op- 
erators of Eqs. (|2.30|) and (|2.31|) . giving results which are free of 0(a) effects at 
all times. In order to increase the signal stability, these correlation functions are 
antisymmetrised with their partners /^^ and /vr- 

Flavour breaking effects have also been monitored by comparing the effective 
pseudoscalar meson mass of two axial correlation functions. The first is the quantity 
defined above, derived from the correlation function composed of a strange (stan- 
dard Wilson) and a down (tmQCD Wilson) quark propagators. The second is the 
effective mass derived from the fully twisted correlator, composed of an up and 
a down (tmQCD Wilson) quark propagator. As the twisted and untwisted quark 
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(L/a)^ X T/a 


2ro 


2^ 


Kg 




^conf 


6.0 


16^ X 48 


0.0931 


1.49 


0.1335 


(0.135169,0.03816) 


402 










0.1338 


(0.135178,0.03152) 


398 










0.1340 


(0.135183,0.02708) 


402 












(U.iooio ( ,U.UZZDi ) 


A nn 
4UU 


6.1 


X 56 


0.0789 


1.89 


0.1343 


(0.1354817,0.0279110) 


100 










0.1345 


(0.1354860,0.0233010) 


99 










0.134/ 


(0.1354896,0.01866 ^8 j 


Izz 


6.2 


243 X 64 


0.0677 


1.63 


0.1346 


(0.1357800,0.0283240) 


200 










0.1347 


(0.1357825,0.0259850) 


201 










0.1349 


(0.1357866,0.0212897) 


214 


6.3 


243 X 72 


0.0587 


1.41 


0.1348 


(0.1358118,0.0246230) 


200 










0.1349 


(0.1358139,0.0222430) 


205 










0.1350 


(0.1358157,0.0198558) 


175 










0.1351 


(0.1358174,0.0174640) 


201 



Table 1: The parameters of the runs at twist angle a = Tr/2. 

flavours have been tuned to be degenerate in mass to O(a^) (cf. subsect. l2.5j) . this is 
a measure of flavour breaking effects. The two effective masses computed from these 
correlation functions differ by ~ 5% at /3 = 6.0, while at /5 = 6.3 they practically 
coincide. Thus, small flavour breaking effects, visible at large lattice cutoffs, disap- 
pear as we move towards the continuum limit. Nevertheless, these effects require 
further and more detailed investigation. 

In order to determine the plateaux of the effective masses we have followed 
the procedure of ref. [31]. We have allowed for a relative excited state contribution 
to the effective mass of at most 0.2% within the plateaux. The plateaux of the 
effective masses [x™"/ro, xjf'^^/ro], which satisfy this criterion, are listed in Table El 
and illustrated in Fig. |2J The values for 3;™™/ro indicate that the K- meson decay 
channel dominates the excited states after roughly 1.25— 1.35 fm from the xq = wall 
source. The pseudoscalar meson mass aMpg is obtained by averaging the aM||(xo) 
values in the plateaux; errors are estimated by the jackknife procedure, omitting one 
measurement from each bin. We present this result in the form rgMps in Table EJ 

Given these considerations, we determine Bk by averaging in the symmetric in- 
terval [x™™/ro, (T — x™™)/ro]. What we compute is the quantity Rb = i?s/ZvA+AV 
(cf. Eqs. ()3. 1913.201 13.211) ). corresponding to the bare Bk-^ We distinguish three 

^Recall that the numerator of Rb is the bare correlation function Fva+av while the denominator 
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sources of uncertainty in our Rb results: 

(i) The results depend on the prescription used for the determination of the nor- 
malisation constants Z\ and Zy, as well as the improvement coefficient Ca- 
As discussed in Appendix IXl two such independent determinations have been 
provided by the ALPHA [38,39] and LANL [40,41] collaborations.'^ We com- 
pute the ratio Rb for both sets of values. The discrepancy between the two 
determinations is a measure of O(a^) systematic effects. We also compute 
Rb with the ALPHA collaboration estimates for Z\ and Zy, but without the 
improving Ca counterterm in the axial current. Any discrepancy between this 
and the previous determinations oi Rb signals the presence of strong 0(a) 
effects. 

(ii) The statistical errors of the correlation functions are computed by a standard 
jackknife error analysis (omitting one measurement from each bin), without 
taking into account the errors of Za, Ca and Zy. 

(iii) The total statistical error is computed by combining the statistical errors of Za, 
Zy (for a given ALPHA or LANL determination) with those of the correlation 
functions, through the standard error propagation procedure. The statistical 
error of Ca is not taken into account, as it is related to an 0(a) correction 
(recall that for the same reason, Csw is also used without an error). 

The results for Rb (and pseudoscalar masses) are collected in Table |2I They 
have all been computed at the time-plateaux shown in the second column, deter- 
mined from the effective mass of Eq. (|4.1|) with ALPHA collaboration values for Za, 
Zy and Ca • The quality of the data is also illustrated in Fig. [2 It is clear that as 
the continuum limit is approached, the O(a^) discrepancies between R^^^^^ and 
^LANL ^gj^(^ iq decrease. This tendency is less marked for R^^^^^ which 
differs from the other two by 0(a) effects. 

The Rb values are extrapolated linearly in (ro-Mpg)^, to the physical point 
{roMxf = 1.5736 (corresponding to = ^[M^o + ^i±]), see Fig. El An in- 
teresting issue concerns the magnitude of the error of these extrapolated values. 
Usually, in quenched simulations performed at a given /3, observables such as the 
ratio Rb are computed on the same configuration ensemble for all values of the 
bare quark masses (hopping parameters). This means that the extrapolation to the 
strange quark mass is performed on data points which are strongly correlated, thus 
reducing the error of the extrapolation. We have instead opted for independent 
Monte Carlo simulations for each set of bare quark masses, which in this respect 

is a properly normalised combination of /a and fv- Also recall that the ratio is improved at tree 
level, and the denominator is 0(a) improved in the chiral limit. 

®The values provided by the LANL collaboration in the more recent ref. [42] only change within 
the quoted errors with respect to the ones we have used in our analysis. 
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™min ™max 
2rn ' 2rn 


roMps 


pALPHA 


pLANL 


jjALPHA; w/o 

Hg 


6.0 


[1.30,3.17] 


2.092(6) 
1.907(7) 
1.780(6) 
1.635(6) 


1.177(9)(8)(12) 
1.139(13)(8)(15) 
1.123(14)(8)(16) 
1.084(13)(8)(15) 


1.089(8)(11)(14) 
1.054(12)(11)(16) 
1.038(12)(11)(16) 
1.001(12)(10)(16) 


1.025(8)(6)(10) 
0.993(12)(6)(13) 
0.977(12)(6)(13) 
0.943(11)(5)(12) 






1.2544 


1.031(19)(23) 


0.952(17)(24) 


0.896(16)(18) 


6.1 


[1.34,3.08] 


1.978(6) 
1.812(7) 
1.647(6) 


1.066(9)(8)(12) 
1.023(10)(7)(12) 
0.981(10)(7)(12) 


1.041(8)(7)(11) 
0.998(10)(6)(12) 
0.955(10)(6)(12) 


0.992(8)(6)(10) 
0.951(9)(6)(11) 
0.911(9)(5)(10) 






1.2544 


0.901(21)(25) 


0.874(20)(25) 


0.835(18)(21) 


6.2 


[1.29,3.05] 


2.079(6) 
1.980(7) 
1.795(7) 


1.023(7)(7)(10) 
0.993(9)(7)(11) 
0.975(8)(7)(11) 


1.017(7)(8)(11) 
0.987(9)(7)(11) 
0.970(8)(7)(11) 


0.980(7)(6)(9) 
0.951(8)(5)(9) 
0.935(8)(5)(9) 






1.2544 


0.902(i3)(31) 


0.900(i3)(32) 


0.8D8(23)(2b) 


6.3 


[1.23,3.00] 


2.050(6) 
1.962(8) 
1.839(10) 
1.722(9) 


0.990(12)(7)(14) 
0.996(15)(7)(17) 
0.959(17)(7)(18) 
0.918(15)(6)(16) 


0.996(12)(5)(13) 
1.003(15)(5)(16) 
0.964(17)(5)(18) 
0.924(15)(5)(16) 


0.962(12)(5)(13) 
0.967(14)(5)(15) 
0.932(16)(5)(17) 
0.892(14)(5)(15) 






1.2544 


0.844(32)(35) 


0.851(32)(34) 


0.820(30)(33) 



Table 2: Results for the pseudoscalar mass and the various ratios from which Bx is 
extracted (twist angle a = 7r/2). The error in r^Mps is statistical. The three errors 
of the Rb ratio are, in order of appearance: (i) due to the statistical fluctuations 
of the correlations; (ii) due to the errors of Z\y, (iii) the total error from the two 
previous ones. The results of the extrapolations to the physical Kaon mass values 
are shown at the bottom of each /3-dataset: the first error Rb is that arising from 
type-(i) errors of the fitted values, while the second is from type-(iii) errors. 

mimic unquenched simulations, at the price of having uncorrelated Rb results with 
larger errors on the extrapolations. 

The justification of the linear fit of i?^, as a function of the squared pseudoscalar 
mass is provided in Appendix O 

4.1 Finite volume effects 

In order to investigate finite volume effects, we have performed a simulation at 
P = 6.0 for a larger spatial volume, with extension L/a = 24, at the lightest mass 
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L3 X T/a^ 


roMps 


pALPHA 


-^conf 


163 X 48 
243 X 48 


1.635(6) 
1.624(5) 


1.084(13)(8)(15) 
1.079(7)(8)(11) 


400 
167 



Table 3: S/^-parameter at fixed effective pseudoscalar mass Mpg for two different 
spatial lattice volumes (twist angle a = vr/2). The various errors on Rb are explained 
in the caption of Table |21 



ro-Mps = 1.635. The results are gathered in Table 01 Finite volume effects are 
estimated in terms of the following ratios: 

Mps(L/a = 16) 000^.3. Rf^^^\L/a = lQ) , _ ^ qq^.^.. .43) 

Mps(L/a = 24) ^ " "-""^^^^ ' i?ALPHA(^/„ ^ 24) ^ " "'""^^l^) • ^4-2) 

Note that all errors are of purely statistical nature. We see that the effective mass 
ratios have a very small statistical deviation from zero, while the Bk ratio is essen- 
tially insensitive to finite volume effects. 



4.2 Non-degenerate masses 



Past and current simulations for the determination of Bk have mostly been carried 
out with degenerate down and strange quarks. The rationale behind this (at least as 
far as quenched simulations are concerned) is to be found in the chiral perturbation 
theory expression for Bk- As discussed in ref. [19-22], the quenched chiral expression 
for Bk is given by 



Bk = B 



1 - {'i + e^)y\ny + hy + cye^ + 5\^- 



2e 



In 



with y = m\/{AT:F)'^ and F ~ 102 MeV (see ref. [43]). The parameter 



md 



(4.3) 



(4.4) 



is a measure of the down-strange degeneracy breaking. The expression (|4.3|) is 
identical in form to the one valid for dynamical quarks, save for the 5— term. This 
quenched artefact vanishes in the degenerate limit e — > 0, but diverges when the 
down quark becomes chiral and e — > 1. 

We have measured Rb with non-degenerate down and strange quarks. Our aim 
is to probe its dependence on quark mass differences, while avoiding the potentially 
dangerous e — > 1 limit. We have performed runs at /3 = 6.0, for two values of 
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e 


Kg 


(Kd,a/i) 


roMps 


pALPHA 


-^conf 


0.00 
0.16 
0.41 


0.1340 
0.1338 
0.1335 


(0.1351830,0.02708) 
(0.1351867,0.02259) 
(0.1351913,0.01598) 


1.780(6) 
1.775(6) 
1.772(7) 


1.123(14)(8)(16) 
1.128(17)(8)(19) 
1.129(14)(8)(16) 


400 
400 
402 



Table 4: i?x-parameter at fixed effective pseudoscalar mass Afpg for three values of 
the degeneracy-breaking parameter e (twist angle a = 7r/2). The various errors on 
Rb are explained in the caption of Table |21 



e ^ 0, tuning the bare mass parameters so that the pseudoscalar mass remains close 
to the value vqMps = 1.780 of the e = simulation. Our results are summarised 
in Table |3 The effect of breaking the mass degeneracy of the valence quarks is 
estimated in terms of the following ratios: 

^^^(^ = °) 1 = 0.003(5), „f^^f(^ = °\- 1 = 0.005(5), (4.5) 



Mps(e = 0.16) Mps(e = 0.41) 



^ALPHA(g^O) ^_ ^B'^™'^(e = 0)__^__^^j^^^^^ 

(4.6) 



0-004(19), „atpha/ - 1 = -0.005(17). 



All errors are statistical. We see that there is no appreciable deviation between e = 
and e 7^ values. It seems that, at least in the region of mass differences explored, 
Bk is not sensitive to the breaking of mass degeneracy. 

Finally, we find it worth to mention in this context a curious by-product of 
our computations that concerns the somewhat surprising appearance of an excep- 
tional configuration in the 7r/2 simulation with non- twisted strange valence quarks, 
corresponding to a rather heavy iT-meson of mx — 650 MeV. This happened at 
P = 6.0 and for a "standard" value of the Wilson hopping parameter which, being 
generally considered "safe" from exceptional configurations, has also been used in 
the simulations of other collaborations. More details are provided in Appendix IdI 



5 Lattice Bk results from tmQCD at a = n/i 

The quenched Bk parameter has also been computed in tmQCD with a twist angle 
a = 7r/4. This formalism, with twisted down and strange quarks in the same flavour 
doublet, is a way to avoid the problem of exceptional configurations, while main- 
taining the property of multiplicative renormalisation of the four-fermion operator. 
Therefore, contrary to the a = 7r/2 simulation, in the present one it is possible to 
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{L/af X T/a 


a 
2ro 


L 

2ro 


{K,an) 




6.00 


24^ X 48 


0.0931 


2.24 


(0.134739,0.010412) 
(0.134795,0.009142) 
(0.134828,0.008397) 


200 


6.10 


24^ X 60 


0.0789 


1.89 


(0.135152,0.00810) 
(0.135190,0.00720) 
(0.135235,0.00615) 


200 


6.20 


32^ X 72 


0.0677 


2.17 


(0.135477,0.007595) 
(0.135539,0.006125) 


73 


6.30 


32^ X 72 


0.0587 


1.88 


(0.135509,0.0076) 
(0.135546,0.0067) 
(0.135584,0.0058) 


76 


6.45 


32^ X 86 


0.0481 


1.54 


(0.135105,0.01459) 
(0.135218,0.01185) 
(0.135293,0.01002) 


105 



Table 5: The parameters of the runs at twist angle a = tt/A. 



simulate at half the value of the physical strange quark mass directly, thus avoid- 
ing the uncertainties introduced by extrapolations from higher masses. The price 
to pay is that larger physical volumes are required at these smaller masses. Our 
simulations have been performed at four values of the lattice spacing in the range 
a ~ 0.05 — 0.09 fm, while lattice spatial extensions are L ~ 1.9 — 2.2 fm. Time 
ranges from xq = to xq = T, with T/L > 2. The two mass parameters (hopping 
parameter k and twisted mass fi for the two degenerate valence flavours) are tuned 
as discussed in sect. Contrary to the 7r/2 case, here we have generated a single 
configuration ensemble per /3, on which measurements at a few {n, afx) values, tuned 
to be close to half the strange quark mass, are performed. The parameters of our 
runs are displayed in Table We point out that at the smallest lattice spacing, 
corresponding to /? = 6.45, APEmille memory limitations made it impossible to 
maintain the physical volume used at the lower (5 values. Consequently, in this case 
we were forced to run at a length of -L ~ 1.5 fm, with larger masses, and extrapolate 
to the pseudoscalar mass at the physical kaon mass. 

The effective pseudoscalar meson mass is computed from the ratio of Eq. 1)4. 1() ; 
again correlation functions are antisymmetrised in time. The excited states analysis, 
carried out as in ref. [44], determines the effective mass plateaux for which the 
relative contribution from the excited states is at most 0.5%. The plateaux of the 



23 



/3 


™min -.max 

^ 2ro 2ro ^ 


roMps 


dALPHA 

Kg 


dL ANL 


j-jALPHA" w/o 


6.0 


[1.30,3.17] 
[1.30,3.17] 
[1.30,3.17] 


1.326(4) 
1.253(4) 
1.207(4) 


1.067(8)(8)(11) 
1.048(8)(7)(11) 
1.036(9)(7)(11) 


0.960(7)(10)(12) 
0.941(7)(9)(11) 
0.930(8)(9)(12) 


0.885(6) (5)(8) 
0.868(7) (5)(9) 
0.857(7) (5)(9) 






1.2544 


1.048(8)(11) 


0.943(7)(12) 


0.869(7)(9) 


6.1 


[1.18,3.55] 
[1.18,3.55] 
[1.26, 3.47] 


1.381(5) 
1.325(5) 
1.257(5) 


0.912(8)(7)(11) 
0.892(9)(6)(11) 
0.865(10)(6)(12) 


0.883(8)(6)(10) 
0.863(9)(6)(11) 
0.836(9)(6)(11) 


0.829 (8)(5)(9) 
0.810 (8)(5)(9) 
0.785 (9)(5)(10) 






1.2544 


0.865(10)(11) 


0.836(9)(11) 


0.784(9)(10) 


6.2 


[1.35,3.52] 
[1.35,3.52] 


1.299(6) 
1.182(6) 


0.882(11)(6)(13) 
0.852(12)(6)(13) 


0.875(11)(6)(13) 
0.845(12)(6)(13) 


0.831 (10)(5)(11) 
0.803 (11)(5)(12) 






1.2544 


0.869(11)(13) 


0.863(10)(12) 


0.818(10)(11) 


6.3 


[1.29,2.94] 
[1.35,2.88] 
[1.47,2.76] 


1.338(9) 
1.259(9) 
1.175(10) 


0.848(13)(6)(14) 
0.831(14)(6)(15) 
0.811(15)(6)(16) 


0.854(13)(5)(14) 
0.835(14)(5)(15) 
0.814(14)(5)(15) 


0.815 (12)(5)(13) 
0.798 (13)(5)(14) 
0.779 (14)(4)(15) 






1.2544 


0.829(13)(14) 


0.834(13)(14) 


0.797(12)(13) 


6.45 


[1.30,2.84] 
[1.25,2.88] 
[1.30,2.84] 


2.054(10) 
1.848(11) 
1.702(11) 


0.937(11)(7)(13) 
0.904(12)(6)(13) 
0.879(13)(6)(14) 


0.949(10)(6)(13) 
0.915(12)(6)(13) 
0.889(13)(6)(14) 


0.913 (10)(5)(11) 
0.881 (11)(5)(12) 
0.856 (13)(5)(14) 






1.2544 


0.821(16)(18) 


0.831(16)(17) 


0.800(16)(17) 



Table 6: Results for the pseudoscalar mass and the various ratios from which Bk is 
extracted (twist angle a = vr/4). The error in roMps is statistical. The three errors 
of the Rb ratio are, in order of appearance: (i) due to the statistical fluctuations 
of the correlations; (ii) due to the errors of Z^y] (hi) the total error from the two 
previous ones. The results of the extrapolations to the physical Kaon mass values 
are shown at the bottom of each /3-dataset: the first error oi Rb is that arising from 
type-(i) errors of the fitted values, while the second is from type-(iii) errors. 

effective masses [x^^'^ / vq, x^^^ / vq] are listed in TableEland illustrated in Fig.|3l Once 
more, they have been obtained with the ALPHA collaboration values for Za, Zy 
and Ca- The quoted values for x™™/ro indicate that the K- meson channel dominates 
the excited states after roughly 1.2 to 1.4 fm from the xq = wall source. The 
pseudoscalar meson mass oMps is obtained by averaging the aMpg(xo) values in the 
plateaux; errors are estimated by the jackknife procedure, omitting one measurement 
from each bin. This result is presented in the form roMps in Table El 

The ratio i?B is computed in the symmetric interval [x^^/ro, (T — XQ'™)/ro]. 
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/3 


L3 X T/a" 


L 

2ro 


T 

2ro 


(k, afj.) 


roMps 


^ALPHA 


N e 


6.0 


lb X 4o 


1.488 




(U. 134739, U. 010412) 


1.316(5) 


1.040(11)(7)(13) 


667 




lb X 4o 


1.488 




(U. 134828, U. 008397) 


1.196(6) 


0.999(13)(7)(15) 


667 




A ri ^ y AO 


2.232 


4.464 


(0.134739, 0.010412) 


1.326(4) 


1.067(8)(8)(11) 


200 




24'' X 48 


2.232 


4.464 


fO 134828 008397) 


1.207(4) 


1 036('9)f7)('l 1) 


200 




32^ X 60 


2.976 


5.580 


(0.134828,0.008397) 


1.201(3) 


1.035(6)(7)(9) 


76 


6.1 


24^ X 60 


1.896 


4.740 


(0.135152,0.008100) 


1.381(5) 


0.912(8)(7)(11) 


200 




24^ X 60 


1.896 


4.740 


(0.135235,0.006150) 


1.257(5) 


0.865(9)(6)(11) 


200 




32^ X 64 


2.528 


5.056 


(0.135152,0.008100) 


1.370(5) 


0.925(7)(7)(10) 


77 




32^ X 64 


2.528 


5.056 


(0.135235,0.006150) 


1.245(5) 


0.876(8)(6)(10) 


77 


6.2 


24^ X 72 


1.620 


4.860 


(0.135477,0.007595) 


1.298(8) 


0.853(16)(6)(17) 


191 




24^ X 72 


1.620 


4.860 


(0.135539,0.006125) 


1.182(9) 


0.813(18)(6)(19) 


191 




32^ X 72 


2.160 


4.860 


(0.135477,0.007595) 


1.299(6) 


0.882(11)(6)(13) 


73 




32^ X 72 


2.160 


4.860 


(0.135539,0.006125) 


1.182(6) 


0.852(12)(6)(13) 


73 



Table 7: Bk parameter (at two values of the effective pseudoscalar mass Mps per 
(3) for different spatial lattice volumes (twist angle a = 7r/4). The various errors on 
Rb are explained in the caption of Table 

The results for Rb are collected in Table El We see once more that as the continuum 
limit is approached, the O(a^) discrepancies between R^^^^^ and R^^^^ tend to 

decrease. This tendency is less marked for R^^^^^ which differs from the 

other two by 0(a) effects. The plateaux from which Bx is extracted are illustrated 
in Fig. ini 

The Rb values are interpolated linearly in (roMps)^, to the physical point 
(ro-Mi^)^ = 1.5736, see Fig. [3 These interpolations are very short and involve data 
which, having been obtained at the same configuration ensemble, are correlated. 
Thus, the error on the interpolated value, compared to the one obtained in the 
a = 7r/2 case by extrapolation from higher pseudoscalar masses, is smaller. For the 
f3 = 6.45 case, we still have correlated data, but as the physical value is reached by 
an extrapolation, the error is relatively large. 



5.1 Finite volume effects 

In order to test finite volume effects, we have performed simulations at /? = 6.0, 6.1, 6.2 
for various spatial volumes and light quark masses. The simulation points, as well 
as the results for the pseudoscalar meson mass and for Rb, are gathered in Table [Tj 
The spatial extent of the lattices considered roughly ranges from 1.5 to 3.0 fm. The 
usual Mpg and R^^^^^ ratios at two different spatial volumes are computed at fixed 
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/3 


a 


L2 
a 


2ro 


L2 
2ro 


(k, afi) 


'Mps(Li) ^ 




ri?ALPHA(^^) 1 
i?^LPHA(^2) 




6.0 


16 
24 


24 
32 


1.488 
2.232 


2.232 
2.976 


(0.134739,0.010412) 
(0.134828,0.008397) 
(0.134828,0.008397) 


-0.008(5) 
-0.009(6) 
0.005(4) 


-0.023(16) 
-0.033(18) 
-0.002(14) 


6.1 


24 


32 


1.896 


2.528 


(0.135152,0.008100) 
(0.135235,0.006150) 


0.008(5) 
0.010(6) 


-0.014(16) 
0.014(17) 


6.2 


24 


32 


1.620 


2.160 


(0.135477,0.007595) 
(0.135539,0.006125) 


-0.001(8) 
0.000(10) 


-0.032(23) 
-0.046(26) 



Table 8: Finite volume effects for pseudoscalar masses Mpg and Bx parameters 
(twist angle a = vr/4). The errors are statistical. 

P and their deviation from unity is presented in Table |B1 

For (3 = 6.0, 6.1 the effective mass ratios display some small finite volume effects, 
at L ~ 1.5 fm to 2 fm, which appear to vanish at L ~ 2.5 fm. For /3 = 6.2 however, 
such effects are also absent already at L ~ 1.5 fm. It has to be noted, however, 
that at this level of precision it is difficult to disentangle finite volume effects from 
e.g. the difference in cutoff effects, or the systematic uncertainties related to excited 
states. 

For Bk we see finite volume effects at L ~ 1.5 fm, which disappear aX L ^ 2 fm. 
This situation is /3-independent. 



6 Bk in the continuum limit 

Having obtained the value of the bare Bk parameter at several lattice spacings, we 
now proceed to the computation of its renormalised counterpart in the continuum 
limit. In ref. [13], the necessary renormalisation constants have been calculated non- 
perturbatively in various Schrodinger functional schemes, as well as the operator step 
scaling functions, which determine its renormalisation group running. What we need 
from ref. [13] is the total renormalisation factor -Hy^_|_^Y.(go)) which connects the 
bare Bk to its RGI value as follows: 

" j™0 ^VA-fAV;s(5o) ^^(So) • (6.1) 

The bare BK{go) is given by Rb at the physical point vqAIk (see Tables |2 and El • 
Which of the three candidates R^^^^^, i^^^^L /j^lphA; w/o ca ^^^^ suitable 
will be discussed below. The quantity ^vA-(-AV-s product of three factors: 

(i) the operator renormalisation constant Zyj^,j^y.^{go,aiJ,aim) at a hadronic scale 
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fJ-min', (ii) the RG-running evolution function from the scale /Umin to a (hopefully) 
perturbative scale /imax = 2^^min; (iii) the RG evolution function from the scale /imax 
to infinity. For several SF renormalisation schemes, the first two factors are known 
non-perturbatively (see ref. [13]), while the third is known in NLO perturbation 
theory (see ref. [45]). Following ref. [13] we write 

^VA+AV;s(5o) = ■Z'vA+AV;s(A'mm) -2^vA+AV;s(50 5 a/^min) , (6.2) 

where .^vA+AV-s(/^min) is the product of the two evolution functions. 

Note that -2^vA+AV s(5'o) depends on the renormalisation scheme s only via cutoff 
effects. In refs. [13] and [45], nine variants of Schrodinger functional frameworks 
have been explored, each giving rise to a distinct renormalisation scheme s. Three 
of these schemes have been deemed "reliable", essentially because they display the 
best control of the systematic uncertainty due to the use of NLO perturbation theory 
in the calculation of -2^va+av s- refer the reader to ref. [13] for details; here we 
just state that, following this work we label these schemes as s = 1,3, 7. 

Given the different origin of the factors contriving to give Bk, the extraction 
of its error is somewhat intricate. The values of -Zy^_(_^Y.g((7o) for the couplings of 
interest are collected in Table IT^ of Appendix where a 1% error is quoted. This 
error is only due to the uncertainties in the determination of the lattice renormal- 
isation constant .^vA+AV-s(50) o^/Umin)- It is combined in quadrature with the error 
for the bare Bxigo) (given in Tables 121 and E]) in order to give the uncertainty for 
Bk, at each j3 value. Our results for Bk (from R^^^^^), with the error obtained 
as described, are gathered in Table IHl Once these results are extrapolated to the 
continuum limit, a further error has to be added in quadrature, due to the evolution 
function Zy^_|^^y.j,(/Xmin); the reason this is only done at the very end is that the 
evolution function is a continuum quantity. As reported in ref. [13], the relative 

error of -^vA+AV;s(^nii° 

) is 2%. 

From the results of Table |H1 we see that Bk is essentially independent of the 
scheme s. In the remaining analysis, we will only consider the s = 1 results. This 
arbitrary choice does not introduce any loss of generality. 

Next we explore the influence of the different choices of axial currents in the de- 
termination of Bk- This entails comparing results extracted from R^^^^^, R^^^^ 
and R^^^^^' see Fig. |H1 Not surprisingly, the 7r/2 data is noisier than that 

of 7r/4. This is due to the fact that the latter dataset has been produced with- 
out extrapolations in the K-meson mass. The comparison between R^^^^^ and 
^ALPHA, w/o CA ii^j-^jg]^ have the same Zj^) shows that there are large CA-related 
cutoff effects at /? = 6.0, which diminish as /3 increases. 

in-between these two cases. Once the f3 = 6.0 data, afflicted by large discretisation 
effects, are conservatively discarded, we are faced with three Bk estimates which 
display reasonable scaling behaviour. At the largest /? values, the three Bk results 
are perfectly compatible. In trying to decide which is the most reliable, we recall 
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13 


6.0 


6.1 


6.2 


6.3 


6.45 


s = 1 
s = 3 
s = 7 


0.911(22) 
0.918(22) 
0.913(22) 


0.812(24) 
0.816(24) 
0.814(24) 


0.828(30) 
0.831(30) 
0.830(30) 


0.789(34) 
0.791(34) 
0.791(34) 


— 


(3 


6.0 


6.1 


6.2 


6.3 


6.45 


s = 1 
s = 3 
s = 7 


0.926(13) 
0.933(14) 
0.929(13) 


0.779(13) 
0.783(13) 
0.781(13) 


0.798(14) 
0.800(14) 
0.799(14) 


0.775(15) 
0.777(15) 
0.777(15) 


0.789(19) 
0.789(19) 
0.790(19) 



Table 9: The value of Bk at fixed /3-coupling, obtained from three distinct 
Schrodinger functional renormalisation schemes s. The quoted error combines in 
quadrature the uncertainty in the computation of the bare -B/^-parameter and that 
of the total renormalisation factor Zyp^_^j^.^{gQ). The upper part of the table refers 
to a = 7r/2 results, while the lower part to a = tt/A. 



that the ALPHA and LANL Za values have been obtained for an axial current 
which contains the Ca counterterm. Thus, the fact that R^^^^' '^^ contains the 
ALPHA Za without the Ca counterterm is somewhat un-natural but not illicit, as 
it only re- introduces 0(a) effects in the axial current. It could even be the source of 
beneficial correlators in the Bk ratio, as the corresponding Ova+av counterterms 
are also missing in the numerator. Nevertheless, the R^^^^' '^^ results would 
have been on a firmer ground, had they been obtained with an independent Za de- 
termination from Ward identities, with the Clover action and without the Ca term. 
Also the R^^^^ dataset is clearly of somewhat inferior quality, as the Za and Ca are 
only given at a few /? values in ref. [41] , while for the other /?'s we had to use simple 
extrapolations (cf. Table [TT|). These general considerations place the R^^^^^ result 
on firmer ground and justify taking the i?^^^^^ dataset (without the /? = 6.0 point) 
as our best candidate for a precision measurement of Bk in the continuum limit. 

A further important element is the choice of extrapolating function. On one 
hand we know that 0(a) effects are present and a linear extrapolation would be 
theoretically justified, which however would amplify the final error. This is especially 
true for our results obtained with a = tt/2, which have rather large errors. On the 
other hand, we see that our data points show an essentially constant behaviour in 
P, supporting a constant fit on empirical grounds, which would give smaller errors. 
The procedure which is theoretically sound, while keeping errors under control, is a 
combined linear fit of the a = 7r/2 and a = 7r/4 results, constraining Bk to a unique 
value in the continuum limit. Following this procedure (see Fig. |^ we arrive at our 
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final result: 



Bk = 0.789 ±0.046. 



(6.3) 



It is nevertheless customary to quote Bk also in the MS scheme at a scale 
2 GeV. The rela ~ 
notation of [13,45]): 



/i = 2 GeV. The relation between Bk and B^^{fi) is known at NLO (we follow the 



Bk 



An 



exp 



'nw 



in 
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libo-lobi 



2bl 



rMS 

Bk 



(6.4) 



The gauge coupling anomalous dimension coefficients bo and bi are universal. The 
LO operator anomalous dimension coefficient ^y^ is also universal, while the NLO one 
7^*" depends on the choice of scheme. For MS one also needs to specify the details of 
the 75 matrix regularisation (i.e. DRED, NDR, HV). The values of these coefficients 
are collected in ref. [45]. Note that in the DRED scheme the renormalised coupling 
we use is the one defined in the MS scheme (rather than the one used in the original 
DRED paper of ref. [46]). Our MS results are 



B 



MS.DRED 
K ' 



2 GeV) 



0.599(36) . 



(6.5) 



B 



MS,NDR 
K 



(2 GeV) = 0.573(34) , 



(6.6) 



^MS,HV^2 GeV) 



0.633(38) . 



(6.7) 



The errors cited above arise from three uncertainties added in quadrature: 

• The 13 K error of Eq. 1)6. 3|) . 

• The uncertainty of , related to the renormalisation group running of g-. 

2.54(8). 



We have used r^A- 



MS 



0.602(48) from [47] to find g^{2 GeV) 



MS" 



• The uncertainty arising from the truncation of the perturbative series, when 
passing from the first, general expression of Eq. (|6.4|) . to the second, perturba- 
tive one. This is estimated by also calculating the perturbative factor through 
numerical integration of the exponent of the general expression (with 7"'" and 
P truncated to the highest available order in PT) . We take the spread of these 
results as our estimate of the 0{g'^) uncertainties. 
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7 Conclusions 



In the present work we have presented a Bk value obtained from first principles, 
without any uncontrolled approximations, except for the quenching of the sea quarks 
and the use of degenerate strange and down valence quarks. The only input from 
experiment has been the setting of the scales from the physical K-meson mass and 
the ro parameter. 

Both matrix element renormalisation and RG-running are non-perturbative. 
NLO perturbation theory has only been used at very high scales, i.e. deep in the 
perturbative regime. The bare Bk has been computed in the lattice regularisation 
with Wilson fermions, using two variants of tmQCD, at several couplings. Thus, a 
well controlled continuum limit extrapolation was possible, in which we have con- 
servatively dropped the result of the coarsest lattice. This was imposed by the large 
cutoff effects which, at small /? values, arise from different choices of the 0(a) im- 
provement term proportional to ca of the axial current. This source of systematic 
error, as well as 0(a) cutoff effects in general, may be removed in future simulations 
performed along the lines proposed in ref. [30]. In any case, our result has been 
obtained with a good control of all sources of uncertainty, except for the quenching 
of the fermion determinant. Most recent quenched results, derived with different 
lattice regularisations, usually with a poorer control of systematic errors, are in 
agreement with our value [4]. It is interesting that this is not the case for the most 
recent measurement of Bk with Wilson fermions [11]. This is discussed in detail in 
Appendix IeI 
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Appendix A Renormalisation constants and mixing coefficients 

In this appendix we collect all the results we have used from other sources, concerning 
bare parameters, renormalisation constants, 0(a) coefficients, matching factors etc. 
Our simulations were produced with the Sheikholeslami-Wohlert action, using 
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the non-perturbative value of the Clover term given in ref. [38]: 



csw - i_ 0.922^2 ' " ^ ^0 < 1 • (A.l) 

For the tuning of the quark masses, as described in subsect. 12. 5| we need the 
values of the critical hopping parameters. These may be found in various sources, 
as shown in Table ITfll 





6.0 


6.1 


6.2 


6.3 


6.4 


6.45 


'^crit 


0.135196 


0.135496 


0.135795 


0.135823 


0.135720 


0.135701 


ref. 


[38] 


[44] 


[38] 


ZeRo Coh. 


[38] 


[44] 



Table 10: The critical hopping parameter Kcrit at several /3-values. 

We also need the non-perturbative expressions for the renormalisation constant 
ZsZp/(ZsZA)ofret. [48], 



Z = (1 + 0.0905149?) ('-°-^'i^«i'g;°-°^;;>;-°-°^373i,g) 



(1 - 0.9678^^ 



, (A.2) 

0<9,^<l 



as well as the axial current normalisation Za, which will be discussed below. At 
present, it suffices to say that we use the ALPHA collaboration estimate of Za for 
the quark mass tuning. Finally, for the mass 0(a) counterterms 6m and 6^ we use the 
non-perturbative and 1-loop perturbative results of refs. [48] and [29] respectively: 



(-0.5 -0.09623,g)(^-";f^°^£t°-°f^^ 



(A.3) 

o<gg<i 



(1 - 0.6905ff, 

= -0.103 Cf gl , (A.4) 

with Cp = {N^ — l)/{2Nc) and Nc = 3 the number of colours. The last expression 
is obtained once the arbitrary value of 6ni has been fixed to ftm = — 1/2 (see ref. [29] 
for a detailed explanation). 

The non-perturbative expressions for the current normalisation constants Za 
and Zy have been computed by the ALPHA collaboration in ref. [39]: 



= l-0.84%,5+ 0.061083 

1 - 0.7332c,2 _ yo - ' v ; 



1 - 0.7663gg + 0.0488gg n ^ 2 ^ «^ 

Zy = K , < Of) < I . (A. 6) 

1 - 0.6369gg - yo - V ; 
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with a precision which is better than 0.6% and 0.4% respectively. The same collab- 
oration finds, for the improvement coefficient of the axial current [38]: 

= -0.007565^ \ ~_ , < < 1 ■ (A.7) 

For the same quantities, the LANL collaboration quote in ref. [41] as their best 
results those reproduced here in Table ITTl In that work, only results at three /3- values 
have been computed. We have interpolated/extrapolated linearly these results (with 
the two independent errors of added in quadrature) in order to obtain estimates 
for all /3-values used in this work. 





6.0 


6.1* 


6.2 


6.3* 


6.4 


6.45* 




0.770(1) 


0.7791(5) 


0.7874(4) 


0.7951(5) 


0.802(1) 


0.8070(9) 




0.802(2)(8) 


0.809(1)(5) 


0.815(2)(5) 


0.818(1)(3) 


0.822(1)(4) 


0.825(1)(5) 


Ca 


-0.038(4) 


-0.036(2) 


-0.033(3) 


-0.033(2) 


-0.032(3) 


-0.031(3) 



Table 11: Current normalisation factors Zy, Za and axial current improvement coef- 
ficient Ca, as reported in ref. [41] at /? = 6.0, 6.2, 6.4. For the other /3- values, indicated 
by an asterisk, we have obtained results by linear interpolation/extrapolation. 

Upon using Z\ and Zy together for the composition of the tmQCD vector and 
axial currents, we account for their correlations by conservatively increasing their 
error to a 2% and 0.5% uncertainty respectively. 

As explained in sect. 13.31 in the computation of Bk the currents we use are 
normalised with Za and Zy and only partially improved by the counterterm Ca- 
The computation of the effective masses oMpg (and eventually the improved decay 
constant Fk) is done using fully improved currents, including the counterterms 
proportional to the quark mass. For the vector current, the counterterm by has 
been computed non-perturbatively in refs. [39,49]: 

l-0.7613gg + 0.0012go^-0.1136g6 n < ^2 < i (A «^ 
1 - 0.91455g ' ^^9o<l, (A.8) 

while the corresponding coefficient 6a of the axial current is given in 1-loop pertur- 
bation theory in ref. [49]: 

6a = 1 + 0.11414 Cf gl . (A.9) 

The twisted mass counterterms for the currents also depend on the arbitrariness of 
6m; setting again b^a — — 1/2? ref. [29] gives at 1-loop in the perturbative expansion: 

6a = 1 + 0.086 Cf gl , (A.IO) 
6v = 1 + 0.074 Cf gl ■ (A.ll) 
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The relation between the lattice spacing and the scale rg is given at ref. [50]: 

ln(a/ro) = -1.6804 - 1.7331(^ - 6) + 0.7849(/? - 6)^ - 0.4428(/3 - 6)^ . (A.12) 

The total renormalisation factor ^va+AV s(5'o)) which connects the bare Bk to its 
RGI value Bk, has been computed non-perturbatively in ref. [13]. In table IT^ we 
collect these results for the couplings of interest. 





6.0 


6.1 


6.2 


6.3 


6.4 


6.45 


^VA+AV;l 


0.884 


0.901 


0.918 


0.935 


0.952 


0.960 


Z+ 


0.890 


0.905 


0.921 


0.937 


0.953 


0.962 


^VA+AV;7 


0.886 


0.903 


0.920 


0.937 


0.954 


0.962 



Table 12: Operator total renormalisation factors ^va_,_av s(5'o) three Schrodinger 
functional schemes (s = 1, 3, 7). The error of these estimates is at most 1%. 



Appendix B Improvement of quenched quark bilinear operators 

In the quenched approximation, the improvement of two-fermion composite opera- 
tors has been carried out in [29] for tmQCD with a degenerate isospin doublet of 
twisted flavours. Here we adapt these results for the cases of interest. The most 
efficient way to do this is by first considering a generalisation of tmQCD, given by 
the action 

Sf = i>{x)[Dw + m + in'y5]^p{x) , (B.l) 

X 

where Dw is the Wilson-Dirac operator, the quark spinor contains in general Nf 
flavours, m and fi are diagonal mass matrices in flavour space and mq the subtracted 
mass matrix. The two actions considered in our work are obtained for Nj = 3 with 
Hiq = diag(0, 0, TTig^s); /X = diag(/i;, — /i;, 0) (for twist angle vr/2) and Nj = 2 with 
mq = diag{mg,mg); fi = diag(/ii, — /i^) (for twist angle vr/4). 

The improvement pattern of a quark bilinear operator Oij = tpiTtpj (with i,j 
denoting distinct flavours and T any Dirac matrix), is established with the aid of 
discrete symmetries of the theory, namely charge conjugation combined with flavour 
exchange rrii mj , <-> Hj , parity combined with sign flips of the twisted masses 
Hi —fii and time reversal combined with sign flips of the twisted masses. Using 



33 



these symmetries we obtain for on-shell improvement 




(B.2) 



Zp[l + bpU 



][Pij - ibpQ 



2 




(B.3) 



2 



(Vr);.,., =Zy[l + 6va^^^^^i^^][y^,,, - i^^a^^l^yl^,,^. + ac^B^T^,] , (B.4) 
(Ar)^,,, = Za[1 + + _ . + acj^p^ , (B.5) 



where 9 is a lattice symmetrised derivative, and T^i, = ^e^^p^Tpfj. The coefficient 
signs are chosen so as to agree with the results of ref. [29] 

Appendix C Linear and chiral log extrapolations 

Here we address the problem of whether it is justified, in the it/2 tmQCD formalism, 
to compute Bk by performing the linear extrapolations of Fig.Elfrom high (roMps)^ 
values to the physical point. An alternative functional behaviour could be that of 
eq. (|4.3|1 . which in the degenerate flavour case becomes 



This is a continuum limit expression. As it is valid close to the chiral limit, there is 
no a priori reason for it to be used in the mass range we are simulating. However, 
its validity can be easily tested on our data. 

We have set three reference values of y, spanning the data points of our 7r/2 
simulations. At fixed /3, we have computed y) for these three reference y- 

values, by linear interpolation (or short extrapolation). We then compute the RGI 
quantity BK{(3-,y) (as explained in sect. IH)) at /? = 6.1,6.2,6.3, for which we are 
confident that scaling has set in. Finally, we extrapolate i?x(/?, y) linearly to the 
continuum limit, thus obtaining three continuum Bxiy) estimates, where Eq. ()C.1|) 
can be applied. 

Next we perform both a linear extrapolation and one based on Eq. (fTTT]) . The 
results of these two extrapolated I3k values are shown in Fig. IIUI where they are 
compared to the continuum Bk estimates, obtained with the same procedure from 
the 7r/4 data, practically computed at the physical mass point. We see that, within 
large errors, the chiral logs are not resolved by our data. Nevertheless, it appears 
that the linear extrapolation performs somewhat better than the chiral log one. 



Bk = B l-3ylny + by . 



(C.l) 
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A last interesting point concerns the comparison of our best Bk estimate (see 
Eq. H6.3() ) to the one obtained in this Appendix, from the Hnear extrapolation: 

= 0.74 ±0.15. (C.2) 

The difference in the extraction of the two results is essentially that the extrapola- 
tions in the mass and the continuum limit have been swapped ^. The fact that the 
two results are compatible means that, for the mass ranges under consideration, the 
order in which these limits are taken is irrelevant. This is clearly not true close to 
the chiral limit. 



Appendix D Exceptional configurations 

The initial motivation behind the introduction of tmQCD was the problem of excep- 
tional configurations [51], which is solved by cutting-off in the IR limit the low- lying 
eigenvalues of the lattice Dirac operator, with the twisted mass fn [12]. The tmQCD 
action we have adopted so far refers to the light flavours, whereas the strange quark 
propagator, being regularised by the standard Wilson fermion action, is not immune 
from exceptional configurations. The common lore [23] is that for the values of Ks of 
Tabled (corresponding to a i^-meson of more than 600 MeV), no exceptional config- 
urations ought to occur. In fact, none of the runs reported above showed any signs 
of such a problem. Nevertheless, evidence of an exceptional configuration showed up 
in an aside simulation, after 363 measurements at (3 = 6.0, a 16^ x 48 lattice volume 
and Ks = 0.1342 (for the record we mention that in the run, = 0.135194 and 
a^d = 0.01051). Thus the problem arose, somewhat surprisingly, at a Kg believed to 
be safe, on empirical grounds. 

In this respect we have considered the effective mass aMpg, as obtained in 
standard (non-twisted) lattice QCD. This is like the quantity computed in Eq. (|4.1j) . 
without the /v correlation functions and with f\ made up of only standard (non- 
twisted) Wilson quarks (i.e. the strange flavours of our formulation). In Fig. ^ 
we plot oAfpg, calculated at the exceptional conflguration as well as at a typical 
configuration of our ensemble. These results are compared to the average oMpg, 
computed on the whole ensemble, save for the exceptional conflguration. This plot 
is a clear warning against simulating even at reputedly safe values of the Wilson 
hopping parameter. We point out that several lattice measurements, besides our 
own, have been performed by other collaborations at this hopping parameter (e.g. 
see ref. [11]). 

^Another important difference is that the value of Eq. 16.31 is the result of a combined fit of the 
7r/2 and 7r/4 data, while that of Eq. I)C.2^ involves only the long extrapolation of the tv/2 data. 
Thus, the errors of the latter estimate are larger. 
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Figure 1: The value of aM^f for: (i) an exceptional configuration; (ii) a typical 
configuration of our ensemble; (iii) the mean value, averaged on the configuration 
ensemble, without the exceptional configuration. 



Appendix E Comparison with the computation of ref. [11] 

In this Appendix we compare our results to the most recent ones obtained with 
Wilson fermions [11], where Bk = 0.96ib0.10 is quoted. There are a few similarities 
and several differences between the two approaches: 

• In both approaches Wilson fermions are used, but in the present work we have 
a twisted term in the lattice action. In ref. [11] the AS" = 2 bare matrix ele- 
ment is computed on periodic lattices, from a four-point correlation function of 
the operator OvA+AV (which is related to the three-point correlation function 
of Ovv+AA through Ward identities). In the present work it is computed from 
a three-point correlation function of the operator OvA+AV (which is related to 
the three-point correlation function of Oyv+AA through the tmQCD formal- 
ism). Schrodinger functional boundary conditions are used, which amounts to 
smeared pseudoscalar sources at the two time boundaries. 

• In both approaches matrix elements are renormalised non-perturbatively. In 
ref. [11] the renormalisation constants are computed in the RI/MOM scheme, 
at scales of a few GeV. The RGI Bk requires also the operator RG-running; 
this is performed in NLO perturbation theory. In the present work, matrix 
elements are renormalised at hadronic scales in a SF scheme. The RG-running 
is performed non-perturbatively, using finite volume techniques, up to a few 
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tens of GeV and subsequently with NLO perturbation theory. 



• In the present work, Bk is extracted directly from suitable ratios of correla- 
tion functions. In ref. [11], Bk is obtained from the slope of the four-fermion 
operator matrix element, as a function of its vacuum saturation estimate 
(8/3)-F|-m|^. This linear behaviour is predicted by lowest order chiral per- 
turbation theory and is also seen in the data produced at pseudoscalar masses 
above the value of the physical X-meson mass. 

It is interesting to try and understand the source of the difference in the I3k 
value quoted by us and by the authors of ref. [11]. First we compare the quan- 
tity -ZvA+AV-i(5o) of ref. [13] to the same quantity, calculated in the the RI/MOM 
scheme; this is done by combining the Zva+av values of ref. [52] , used in the Bk com- 
putation of ref. [11], with the NLO running quoted in refs. [2,3]. The two estimates 
do not depend on the renormalisation scheme in which they have been computed. 
They only differ by (i) discretisation effects; (ii) the reliability (convergence) of the 
NLO perturbative running, calculated in two different schemes and applied at two 
different scales. In Fig. ^2 we show the quantity 

^ _ ^VA+AV;l ^VA+AV;RI/MOM 
^VA+AV;1 

as a function of the lattice spacing. It is clear that a barely significant difference at 
the largest lattice spacing, quickly disappears as the continuum limit is approached. 
Thus it seems that renormalisation is not responsible for the different Bk results 
quoted by the two groups. 

In order to have as direct a comparison as possible, we compute in the a = it/2 
case the quantity 

j^r ^ ^-^VA+Avi^^A+AV ^-^ 2) 

2[/p(xo) - i^fsixoWUT - xo) - i^niT - xo)] 

Z-p Z:p 

The computation is performed at each value of the inverse coupling (3 and mass 
parameters Kg and {n^^a^) of Tabled The normalisation factor Zo,/Z-p = {ZZa)~^ 
is calculated from the Z and Za ALPHA collaboration estimates quoted in Ap- 
pendix ^ 

A glance at Eqs. (12.1112.121 IB. 2|B~3|) shows that the above ratio tends, at large 
time separations, to the quantity 

^, ZvA+Av(i^°|QvA+Av|i^°) 

mPs,\K')\' ■ ^ ^ 

Finally, for each set of coupling and mass parameter values, we construct the quantity 

-^RGI = ^VA+AV;l(5'o) R' (E.4) 
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in exactly the same way we have computed Bk of eq. (|6.1j) . This quantity may be 
directly compared to the ratio of ref. [11] (see Table 1 of that work; the 

same quantity is sometimes referred to as R''^'' by the authors).^ 
We also compute the quantity 

X' = ^ ^■^^■^^(^"^ ~ ^■^v/v(^o)][^A/p(r - xq) - iZyf^{T - xq)] ^^^^ 
3 [/p(:ro) - ^^f,{xo)mT - xo) - z^fi{T - xo)] 

Z-p Z-p 

This ratio tends, at large time separations, to the quantity 

3K0|P,,|K0)P' ^^•'^^ 

which is the same quantity as X of ref. [11]. For both R'^qi and X' we have used 
ALPHA collaboration renormalisation constants. 

At fixed gauge coupling /3, the quantity R^qi is fit linearly as a function of X'; 
the slope is the Bk estimate according to lowest order chiral perturbation theory. In 
Fig. ^2 we compare the results of ref. [11] to our results (obtained as described here 
from the ratio -Rj^qj) and find them compatible. We also see that our preferred Bk 
estimates (obtained on the same gauge configurations, but extracted from the ratio 
Rb) give a consistently lower Bk estimate. This shows that the discrepancy of the 
two results is not due to different regularisation and renormalisation systematics, 
but due to the different method of extraction of Bk- 

The errors we quote for the results of the present method are larger than those 
of ref. [11]. This is due to the fact that the latter results have been obtained, at 
fixed P and for several values of the hopping parameters, on the same configuration 
ensemble, while we have generated a separate ensemble at each k. Thus, the extrap- 
olations in the pseudoscalar mass in ref. [11] benefit from the correlations between 
measurements. On the other hand, for both datasets, Bk has been computed from 
the slope of R'^qi as a function of X'. This fitting procedure enlarges the errors, as 
can be seen by comparing them to those of our main method (based upon computing 
-Rb at several k- values and extrapolating it to the physical Kaon mass). 
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Figure 2: Plateaux for the extraction of Mpg at a = 7^/2. The time-range and value 
of each plateau is indicated by a straight line segment. 
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Figure 3: Plateaux for the extraction of Mpg at a = 7r/4. The time- 
of each plateau is indicated by a straight line segment. 
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Figure 4: Plateaux for the extraction of R^^^^^ at a = 7r/2. The time-range and 
value of each plateau is indicated by a straight line segment. 
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Figure 5: Plateaux for the extraction of R^^^^^ at a = 7r/4. The time-range and 
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Figure 6: Linear extrapolation of Rb to the physical kaon mass at a = tt/2. 
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Figure 7: Linear inter /extrapolation of Rb to the physical kaon mass at a = tt/A. 
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Figure 9: Continuum limit extrapolation of Bk, obtained from a combined linear 
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